Cavity opto-electromechanical system combining strong electrical actuation with ultrasensitive 

transduction 
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A cavity opto-electromechanical system is reported which combines the ultrasensitive transduction of cavity 
optomechanical systems with the electrical actuation of nanoelectromechanical systems. Ultrasensitive mechan- 
ical transduction is achieved via opto-mechanical coupling. Electrical gradient forces as large as 0.40 /uN are 
realized, facilitating strong actuation with ultralow dissipation. A scanning probe microscope is implemented, 
capable of characterizing the mechanical modes. The integration of electrical actuation into optomechanical 
devices is an enabling step towards the regime of quantum nonlinear dynamics, and provides new capabilities 
for quantum control of mechanical motion. 



I. INTRODUCTION 

Mechanical oscillators offer the potential to test quan- 
tum mechanics at the macroscopic scale |[l|-[3], as well as 
significant technological advances including applications in 
metrology 1 7] and quantum information systemsisl [2]. Two 
key requirements to enter this quantum regime are the ca- 
pacities to efficiently quantum control the oscillator, and to 
achieve zero-point motion limited transduction sensitivity. 
Strong electrical actuation of mechanical oscillators, an en- 
abling step for quantum control, has been achieved in nano- 
electromechanical systems (NEMS) |[1Q| - [15|] . Eurthermore, 
transduction sensitivities near the quantum limit have been 
reached with cavity optomechanical systems (CQMS) |[l6l 
nm . We have recently combined these two techniques into 
the first cavity opto-electromechanical system (C0EMS) |[l8ll : 
achieving gradient forces an order of magnitude stronger than 
any cryogenic radiation pressure based microcavity actuation 
technique to date, with intracavity optical enhancement en- 
abling transduction at levels close to the mechanical zero- 
point motion iIitI] . In this paper we present further analysis 
of this system including new theoretical and experimental re- 
sults confirming gradient forces as the mechanism responsible 
for mechanical actuation, an explicit derivation of the peak 
to peak gradient force and new results imaging mechanical 
modes with COEMS. 

The capacity to strongly electrically actuate COMS re- 
ported here has several important ramifications. Electrical 
feedback of the optomechanical transduction signal allows 
immediate control of the properties of the mechanical oscil- 
lator; facilitating, for example, active cooling or heating 0- 
[20], electro-optic spring effects! 21], and phonon lasing|22]. 
Strong electrical driving of the mechanical motion provides 
access to mechanical nonlinear dynamics in both classical and 
ultimately quantum regimes; with nanofabrication techniques 
providing the means to engineer the nonlinear properties. 
This nonlinear mechanical regime is usually achieved with 
NEMS. Eor example the coherent response of mechanical os- 
cillations allow signal enhancement and noise reduction 11231] . 
ultrasensitive mass detection iH U^, sub-Heisenberg limit 
metrology 1 3], and mechanical quantum state engineering IJ]. 
The superior transduction sensitivity afforded by COMS 



presents an enabling step towards observation of quantum 
nonlinear dynamics in these systems. Eurthermore, the inte- 
gration of electrical actuation into COMS advances the possi- 
bility of incorporating ultracold superconducting circuits, and 
unifying cavity opto-mechanics with the rapidly developing 
field of circuit quantum electrodynamics 11261] . 

COMS consist of a high quality optical cavity strongly 
coupled to a mechanical oscillator. Motion of the mechan- 
ical oscillator alters the cavity's optical path length, and as 
a result this motion is transferred to the out-coupled optical 
field. Information extracted from measurements of this field 
enables both direct characterization of the mechanical oscil- 
lator, and, combined with mechanical actuation, control of its 
behaviourlUsl |20] . To date, the actuation has been achieved 
using radiation pressure from an amplitude modulated opti- 
cal field (for example see Refs. |[l9ll20,l22l])- However, ra- 
diation pressure is an inherently weak force. Eurthermore, 
optical heating via optical absorption of intracavity power 
severely constrains the possibility of operation in the quan- 
tum regime ! 171 . By contrast, electrical actuation of NEMS 
is inherently strong, much less prone to heating effects, and 
comparatively straightforward to implement |[lll - [l4l \2M . 

NEMS are submicron structures that use actuators to con- 
vert electrical energy to mechanical energy by exciting a 
resonant mode structure. These structures have fundamen- 
tal resonance frequencies in the radio frequency (RE) to mi- 
crowave frequency domain, mechanical quality factors ,Qm, 
in the hundreds of thousands and masses in the femtogram 
range 1123]. External actuation mechanisms that use dieletric 
gradient forces reduce constraints on material choice and sys- 
tem geometry, and also avoid dissipative actuation losses lllSll . 
These dieletric gradient forces are the result of a non-uniform 
applied electric field that polarizes the dielectric so that it ex- 
periences a dipole force in the direction of increasing electric 
field. NEMS, however, have relatively poor transduction sen- 
sitivity on the order of 10"^^ m Hz~^^^ II3Q| - [33|] in comparison 
to the COMS mentioned above which achieve transduction 
sensitivities better than 10"^^ m Hz"i/^|[i7|]. The COEMS pre- 
sented in Ref. lflSn integrates dieletric gradient force actuation 
with silica microtoroiods on a silicon chip which have both 
high quality mechanical resonances and ultra high quality op- 
tical resonances. This simultaneously allows gradient forces 



as large as 0.4//N, and ultrasensitive transduction at the level 
of 1.5x10-1^ mHz-i/2. 

To demonstrate the broader relevance of this work, we im- 
plement COEMS-based scanning probe microscopy of the mi- 
crotoroid mechanical mode structure. This is a natural ex- 
tension of the ability to both precisely actuate and measure 
the displacement of the resonator. Microscopy based on the 
measurement of a displacement of some spring structure in 
response to a force is now common-place. Atomic force mi- 
croscopy (AFM) ll34ll is able to produce three dimensional im- 
ages with atomic scale resolution of both conductors and insu- 
lators. The force on the probe may be due to a range of prox- 
imity forces including van der Walls, capillary, electrostatic, 
and Casimir forces; and is kept small and constant with feed- 
back mechanisms 11351 [3^. In another example, magnetic res- 
onance force microscopy is based on the detection of the mag- 
netic force between a ferromagnetic tip and the electron spins 
in a sample and can detect single electron spins 11371] . Of more 
direct relevance to the work reported here, the spatial profile of 
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FIG. 1 : (color online) A Schematic of gradient force actuation ex- 
periment. FPC: fibre polarization controller. Inset: electric field dis- 
tribution between sharp and flat electrodes. B The method of im- 
ages used to describe the electric field between the sharp and flat 
electrodes. Left to right: mechanical and electrical structure; struc- 
ture showing induced charges in sharp and flat electrodes, and mi- 
crotoroid polarization; structure showing electric field lines modeled 
using the method of images where the flat electrode can be modeled 
as an equivalent but opposite sign charge to that on the electrode tip 
located equidistant beneath the flat electrode itself. 



acoustic modes in high finesse Fabry-Perot cavities have been 
imaged where the actuation force has been achieved via both 
optical radiation pressure[38] or electrostatic actuation 11391] . 
In this paper we extend these techniques to measure a two 
dimensional spatial profile in integrated microtoroid cavities. 
Although currently limited in resolution due to the size of the 
region which our actuation force is applied to; this technique 
allows confirmation of the mechanical mode shapes of the mi- 
crotoroid and gives a spatial profile that is in agreement with 
finite element modeling (COMSOL Multiphysics). With the 
implementation of feedback similar to that used in AFM to 
control the probe position sub-micron resolution should be 
possible. 



II. EXPERIMENT 

Figure [TK shows a schematic of our experimental setup. 
The microtoroids were fabricated on a silicon wafer with a 2 
yum layer of wet grown silica. The silica was patterned us- 
ing standard photolithography techniques, followed by etch- 
ing in hydrofluoric acid to transfer a pattern of 100 yum diam- 
eter disks into the silica. Xenon difluoride was then used to 
isotropically etch 25 yum of silicon enabling the undercut sil- 
ica disk to be reflowed with several watts of focused CO2 laser 
power 1140] . The microtoroids then typically had major and mi- 
nor diameters of 60 yum and 6 yum, respectively, and an under- 
cut of 10 yum between the edge of the toroid and the pedestal. 
An external cavity diode laser (New Focus TLB-6320) that 
could be tuned over several microtoroid free spectral ranges 
provided an in-fiber optical probe field at 980 nm. The probe 
field was evanescently coupled into a microtoroid using a sin- 
gle mode tapered optical fiber held within several micrometers 
of the microtoroid surface with a three axis piezo-driven flex- 
ure stage. 

We characterize the optical quality factor ,2^, by fit- 
ting the observed transmission spectra to a coupled cavity 
modelllll |42]. The intrinsic quality factor was found to be 
2o = 3x10^ enabling thermal locking of the microtoroid to the 
probe laser. This thermal locking is a result of the temperature 
dependent refractive index of silica which with optical heating 
shifts the microtoroids resonant frequency to remain locked to 
the probe laser ll43l] . All experiments were carried out at room 
temperature {T = 300 K) and atmospheric pressure. 

To achieve dielectric gradient force actuation a non-uniform 
electric field was produced with a RF voltage from a net- 
work analyzer (Hewlett Packard 3577B) and applied to a sharp 
stainless steel electrode. This sharp electrode had a 2 yum 
tip diameter and was positioned 15 yum vertically above the 
center of the microtoroid, with the height chosen to balance 
the benefit of increased gradient forces against the risk of 
chemical contamination and structural damage due to phys- 
ical contact. Stepper motor drives (Thorlabs DRVOOl) en- 
abled both accurate positioning and fast scanning of the sharp 
electrode. The silicon wafer with integrated microtoroids was 
placed between this electrode and a large flat electrode that 
was grounded to a voltage supply but could be oflTset for par- 
ticular measurements. A schematic of the electric fields pro- 



duced by our electrode arrangement can be seen in Fig. [T^ 
where the method of images, typically used to describe bound- 
ary conditions in electrostatics, is used to predict the form of 
the electric field lines. Figure [T^ shows that a voltage applied 
across the sharp and flat electrodes induces charge separation 
and the electric field through the microtoroid is well approxi- 
mated by that of a point charge located at the tip of the sharp 
electrode. This was confirmed with finite element modeling. 
As will be shown later, our experiments suggest that the in- 
trinsic polarization of the silica microtoroid material due to 
static electric fields dominates the induced polarization due to 
our RF actuation. This induced polarization is most probably 
caused by charges trapped in the silica. The inhomogeneous 
electric field acts on the polarized silica creating the actuation 
force to stimulate the mechanical modes of the microtoroid. 



III. RESULTS 
A. Ultrasensitive mechanical transduction 

Detuning the probe field to the full-width-half-maximum 
(FWHM) of the optical mode allowed transduction of the me- 
chanical motion; with mechanical motion induced shifts in the 
microtoroids optical path length being transferred to the am- 
plitude of the out-coupled optical field and detected directly 
with a silicon photodiode (New Focus 1801)|Hl- The trans- 
duction noise spectral density obtained by spectral analysis of 
the output photocurrent, without application of an RF drive 
voltage, is shown in Fig. |2] for 1 mW of incident probe power. 
By applying a known reference modulation we were able to 
calibrate the absolute amplitude of the displacement from the 
optomechanical transduction signal [14511. The three spectral 
peaks in Fig. 2 are due to the Brownian motion of the mi- 
crotoroid mechanical modes. The peaks were identified with 
finite element modeling as the lowest order flexural mode and 
the two lowest order crown modes, where the frequencies of 
the modes and the model agree to within 5%. The flat spectral 
noise background was due to laser phase noise, as confirmed 
by an observed square root dependence of its amplitude on 
probe power. This noise placed an ultimate limit on the trans- 
duction sensitivity. 

The quantum behavior of a mechanical oscillator can only 
be observed if the transduction sensitivity is sufficient to re- 
solve the oscillators zero point motion. Brownian motion can 
be modeled as the mechanical response due to a Langevin 
force Ft that describes the coupling of the oscillator to the 
thermal bath with thermal energy |^r7j20]. The spectral den- 
sity then has a Lorentzian lineshape 11461] 

sf{co)W}lz~'] = \)(j{co)\^SF^=2kBTYjmj\)(j{co)\\ (1) 
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valid when the spectral analysis is at frequencies close to me- 
chanical resonance. In our case the total measured spectral 
density is the sum of the spectral density of the Brownian 
motion of each mechanical mode, and a spectral background 
noise ^a^ due to transduction noise which is approximately flat 
in frequency 
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Using Eqs. (lll2l3l) it was possible to accurately fit the observed 
mechanical spectral density as shown in Fig|2l with ntj, Tj, 
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and ^A^ as free parameters. From this fit we find S^ 

2Q-18 j^ Hz"^/^, and the effective masses and damping rates 

for each of the three modes shown in Table U 





TABLE I: Mechanical mode parameters 
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The peak of the spectral density of a mechanical mode due 
to zero-point motion is given by ll45n 
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which gives values in the range of 10"^^ m Hz"^/^ for the 
modes considered here, as shown in Table U Thus our trans- 
duction sensitivity is approximately two orders of magnitude 
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Here, the parameter j identifies each mechanical mode with 
7= 1 being the lowest order crown mode, j=2 the lowest order 
radial ffexural mode and j=3 the second order crown mode; 
and Mj and ojmj are the effective mass and resonance fre- 
quency of mode j respectively. The dissipation rate Yj is the 
FWHM of the power spectrum and the mechanical suscepti- 



FIG. 2: (color online) Square-root Brownian motion displacement 
spectral density. Black curve: measured spectra; solid red (grey) 
curve: theoretical model; dashed red (grey) line: transduction sen- 
sitivity limit. Insets: (from left to right) finite element models of 
j - 1, 2, 3 modes respectively. 
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FIG. 3: (color online) Gradient force actuation of a COEMS. A 
Square root transduction spectral density. Black curve: measured 
spectra; red curve: theoretical model including interference between 
the mechanical modes; dash-dotted line: transduction sensitivity 
limit. B Square root peak transduction spectral density as a function 
of RF drive amplitude. Green triangles: 7 = 1, red circles: j = 2, 
black squares: 7 = 3. 



away from the zero-point motion as shown in the last column 
of Table II The best transduction sensitivity achieved is just 
over thirty times the zero point motion. 

Recently a number of techniques have been developed to 
improve micro toroid based COMS. These include homodyne 
and polarization spectroscopic techniques that allow the de- 
tection of the phase of the light with quantum limited sensitiv- 
ity, and have recently been applied to improve the transduction 
sensitivity for motion transduction in WGM resonators ll45ll . 
Reduction of the mechanical mode damping rate has been 
achieved through microfabrication ll47ll . Cryogenic cooling 
limits noise due to temperature induced fluctuations in the re- 
fractive index and placing the microtoroid in vacuum reduces 
the mechanical damping due to air pressure 11481] . Mechanical 
Q factors of 80 000 have been reported with these techniques 
while retaining ultrahigh optical Q. The COEMS discussed in 
this paper is compatible with these techniques and their imple- 
mentation should enable zero-point motion transduction sen- 
sitivity. 



B. Driving mechanical modes 

Electrical gradient force actuation of COEMS ofl'ers the 
possibility of broadband excitation, is compatible with inte- 
grated circuits, and is scalable to large architectures. The gra- 
dient force actuation of our COEMS was characterized with 
a network analyzer. The output voltage from the network an- 
alyzer was supplied to the sharp electrode and the frequency 
response was used to monitor the transduction signal. The 
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FIG. 4: (color online) Characterization of the COEMS gradient 
force. A Square root peak transduction spectral density for fixed 
RF drive amplitude as a function of applied DC voltage. B Power 
spectrum peak as a function of probe height z. 



observed mechanical frequency response is shown in Fig.[3j\ 
for a 3 Vrms network analyzer output voltage. A large increase 
in the mechanical oscillation was observed as the RF signal 
was scanned over the mechanical resonance frequency. The 
inverted part of the lowest order crown mode is due to inter- 
ference between the two lowest order modes as a result of the 
coherent driving field. At 3 Vrms the maximum oscillation 
amplitude observed was S J^^^^^ = 2.4 x 10"^"^ m Hz"^/^ at the 
peak of the second order crown mode (7=3). The peak to peak 
gradient force can be derived from this measurement in the 
following way. The dynamics of a damped driven harmonic 
oscillator can be modeled by the equation of motion 
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where x is the displacement of the oscillator from equilibrium, 
k = o/:mj is the spring constant, and fo is the actuation force. 
Taking the Fourier transform and rearranging we find 
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The root-mean-square force is then 
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where we have used the fact that {x(ojj)^) = Sx(coj)Trbw- The 
peak-to-peak force is simply related to the rms force by Fpp = 
2 ^/2FRMs so that 
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For the second order crown mode measured above, we find a 
peak to peak gradient force of 0.40 yuN, surpassing all cryo- 
genically compatible COMS to date by more than an order 
of magnitude. The linear relationship between the applied RF 
voltage and the peak oscillation amplitude of each mechanical 
mode is shown in Fig. [3^ where we see the displacement is 
proportional to the peak drive amplitude (Eq.[8]). 

Gradient force actuation relies on the presence of a net 
dipole moment. The dipole moment, p, of the silica in our 
microtoroid has an intrinsic component due to charges trapped 
in the silica and an induced component due to the applied DC 
field, that is 



P — /^intrinsic "I" /^induced- 



(9) 



The electric field of a dipole along the vertical axis z is given 
by 
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where p^ is the z component of the total dipole moment p. As 
the electric field has a z^ dependance, the force exerted by the 
dipole (the microtoroid) on a point charge (the sharp probe) is 
expected to also go as z^. Newton's third law tells us that the 
force exerted by the point charge on the dipole must be the 
same, and therefore also go as z~^ • Furthermore, since ^induced 
is proportional to the applied DC field, the applied force is 
expected to vary linearly with DC voltage. 

The response of our COEMS to applied DC voltage, as 
predicted above, allowed confirmation that the applied forces 
were indeed due to gradient force actuation. The peak oscil- 
lation amplitude was found to vary linearly with DC voltage 
(Fig.IlK) as expected. Furthermore measurements with vary- 
ing electrode height z confirmed this demonstrating that the 
oscillation amplitude decays as z^ as shown in Fig. |4^ as 
predicted. Note that by selecting the DC voltage such that the 
induced polarization has equal magnitude but opposite sign to 
the intrinsic component it is possible to zero the z-component 
of the net dipole moment and hence switch off' the RF fre- 
quency response. By selecting the correct DC voltage can the 
gradient force can be turned on or off' and hence act as a switch 
in photonic circuits ll49ll . 

There was no observable degradation in the mechanical 
quality factor of the microtoroid caused by the sharp electrode 
or the RF drive voltage indicating that the actuation mecha- 
nism is low dissipation. To quantify this eff'ect we measure 
the mechanical dissipation rate at a range of drive voltages 
for a typical microtoroid mechanical mode. Figure [5]\ shows 
the dissipation rate for peak spectral densities varying from 
8 to 120 times the peak of the background Brownian noise. 
Lorentzian fits allowed the dissipation rate to be extracted in 
each case, with the results plotted in Fig. [5^. As can be seen, 
within experimental uncertainties the dissipation rate is unef- 
fected by the drive voltage amplitude; with the line of best fit 
to the dissipation rate versus spectral peak amplitude having a 
slope of 40+500 mHz/dB 

For quantum behavior to be observed in a COMS requires 
temperatures T <^ Hojmj/kB. For a mechanical resonance of 
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FIG. 5: (color online) Characterization of the dissipation of the CO- 
EMS actuation. A Mechanical transduction spectra with varying 
drive amplitude. The highlighted strip indicates the FWHM. B Dis- 
sipation rate, F, for each drive amplitude including a line of best fit. 



say 5.5 MHz this requires a very low temperature of 0.26 mK. 
The ultra-high Q cavities required for radiation pressure ac- 
tuation experience substantial heating due to intracavity ab- 
sorption constraining the possibility of reaching the quantum 
regime. It has been shown that 1 W of intracavity power 
causes a lOK temperature increase in microtoroids similar to 
those used in this work lflTI] . With radiation pressure alone an 
intracavity power of P^ cF^ms /tt = 36 W would be required to 
produce the maximum oscillation amplitude demonstrated in 
this work ll5Qll . These power levels therefore clearly cause suf- 
ficient heating to be incompatible with the quantum regime. 
The capacity of our COEMS to achieve transduction sensitiv- 
ity close to the standard quantum limit and and actuation with- 
out appreciable heating should facilitate future quantum com- 
munication and quantum nonlinear mechanical experiments. 



C. Microscopy of mechanical mode structure 

Displacements in a microtoroid can result from either the 
propagation of mechanically resonant modes that deform the 
microresonator structure, as discussed above, or from dis- 



placement of the entire resonator due to environmental fac- 
tors. Movements of the entire structure are not the focus of 
this paper as they are strongly dependent on how well the 
microtoroid is isolated from environmental perturbations and 
usually have low resonance frequencies that can be easily dis- 
tinguished from the internal resonant modes Il38ll . As the qual- 
ity of the finite element modeling used to predict both frequen- 
cies and modeshapes of resonant modes depends heavily on 
the initial mesh chosen it is important to validate the simula- 
tions with experimental results. Here we demonstrate a simple 
technique to directly measure the driven mechanical modes as 
previous techniques for measuring mechanical modes Il38i 13911 
are not directly applicable to these whispering gallery mode 
resonators. 

Electrostatic actuation has been previously applied II39[1 as 
the actuation force to form a cross sectional image of the me- 
chanical modes of Fabry-Perot resonators. Here we apply the 
electric gradient force to integrated microtoroid cavities, by 
automatically scanning the sharp electrode across the micro- 
toroid surface while strongly electrically driving the system at 
a mechanical resonance frequency and using optomechanical 
transduction to observe the mechanical response. This allows 
a two dimensional picture of the resonant mode to be gener- 
ated. 

The sharp electrode was positioned 15 jim above the mi- 
crotoroid surface and scanned laterally in the XY plane using 
computer controlled stepper motor drives. A fixed-frequency 
drive voltage was applied with frequency chosen to coincide 
with the peak oscillation frequency of the microtoroid radial 
flexural mode, and the driven peak oscillation amplitude was 
observed throughout the scan using opto-mechanical trans- 
duction. Consistent with this choice of mode, a radially sym- 
metric profile was found, as shown in Fig. [6]\. In contrast, 
crown modes, whose mechanical motion exhibits sinusoidal 
oscillations around the circumference of the torus ll45ll , are ex- 
pected to have a minimum at the microtoroid center, since the 
forces on out-of-phase nodes cancel; and this was indeed ob- 
served, as shown in Fig. [6^. The spatial resolution was lim- 
ited by the non-point-like nature of the electric gradient force; 
with the observed structure being a convolution of the actual 
mechanical mode structure and the spatial distribution of the 
gradient force. This limitation could be overcome in future by 
positioning the sharp electrode closer to the microtoroid sur- 
face, using electronic control loops similar to those in atomic 
force microscopy 1I34I1 to avoid unintentional contact. 



IV. CONCLUSIONS 

The COEMS architecture described in this paper combines 
the ultrasensitive transduction capabilities of COMS with the 
gradient force actuation capabilities of NEMS. Forces as large 
as 0.40 jjN were applied to the mechanical resonator due to 



the inhomogenous electric field from a sharp electrical probe 
tip. These actuation forces do not result in significant heating 
and have significantly higher magnitude than what has been 
achieved with cryogenically compatible radiation pressure ac- 
tuation. This COEMS system has a transduction sensitivity of 
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FIG. 6: (color online) Scanning probe microscopy of mechanical 
mode structure. A Square root peak transduction spectral density 
for a radial flexural mode as a function of the lateral position of the 
sharp electrode. B Cross-sections of the square root peak transduc- 
tion spectral density for a radial flexural mode (dark blue trace) and 
a second order crown mode (light red trace). 



1.5 X 10"^^ m Hz"^/^ on the order of thirty times the mechani- 
cal zero-point motion. Lateral scanning of the drive electrode 
enabled microscopy of the mechanical mode structure, pro- 
viding a tool to facilitate quantitative understanding of the 
mechanical properties of optomechanical systems. Measure- 
ments of the change in transduction signal with the probe 
height and applied DC voltage were used to show that the 
intrinsic polarization of the silica was the origin of the res- 
onators mechanical response to the gradient field. By apply- 
ing a DC bias voltage it was possible to turn the mechanical 
resonators response on and off', giving the COEMS a photonic 
switch capability. 

This research represents important progress in our capac- 
ity to control cavity optomechanical systems at the quantum 
level, allowing, for example, feedback cooling and electro- 
optic spring eff'ects; and presenting an enabling step to- 
wards the new regime of quantum nonlinear mechanics, where 
strong mechanical driving of a ground- state cooled mechani- 
cal oscillator allows exploration of quantum nonlinear dynam- 
ics and nanofabrication techniques provide the capacity to en- 
gineer the nonlinear properties. 

This research was funded by the Australian Research Coun- 
cil grant DP0987146, and performed in part at the Australian 
National Fabrication Facility. We gratefully acknowledge 
helpful advice from T. Kippenberg, G. Milburn, and T. Stace; 
and experimental support from A. Rakic and K. Bertling. 



[1] K. C. Schwab and M. L. Roukes, Physics Today 58, 36 (2005). 



[2] D. Rugar and R Grutter Phys. Rev. Lett. 67 699 (1991). 



[3] M. J. Woolley, G. J. Milburn and C. M. Caves, New J. Phys. 10, 

125018 (2008). 
[4] I. Katz, A. Retzker, R. Straub and R. Lifshitz, Phys. Rev. Lett. 

99 040404 (2007). 
[5] W. Marshall, C. Simon, R. Penrose and D. Bouwmeester, Phys. 

Rev. Lett. 91 130401 (2003). 
[6] A. D. OConnell, M. Hofheinz, M. Ansmann, R. C. Bialczak, M. 

Lenander, E. Lucero, M. Neeley, D. Sank, H. Wang, M. Weides, 

J. Wenner, J. M. Martinis and A. N. Cleland Nature 464 697- 

703 (2010) 
[7] C. L. Degen, M. Poggio, H. J. Mamin and D. Rugar, Phys. Rev. 

Lett. 100, 137601 (2008). 
[8] S. Mancini, D. Vitah , P Tombesi, Phys. Rev. Lett. 90 137901 

(2003); 
[9] P. Rabl, S. J. Kolkowitz, F. H. Koppens, J. G. E. Harris, P. 

Zoller, M. D. Lukin, arXiv:0908.0316 (2009). 
[10] J. D. Teufel, J. W. Harlow, C. A. Regal, &K. W. Lehnert Phys. 

Rev. Lett. 101, 197203 (2008). 
[11] T. Rocheleau, T. Ndukum, C. Macklin, J. B . Hertzberg, A. A. 

Clerk and K. C. Schwab iarXiv:0907.33 131 /1 [quant-ph] July 

(2009). 
[12] M. Poggio, C. L. Degen, H. J. Mamin, & D. Rugar Phys. Rev. 

Lett. 99 017201 (2007). 
[13] D. G. Blair, E. N. Ivanov, M. E. Tobar, P J. Turner, R van Kann 

& I. S. Heng Phys. Rev. Lett. 74 1908-1911 (1995). 
[14] J. B. Hertzberg, T. Rocheleau, T. Ndukum, M. Savva, A. A. 

Clerk and K. C. Schwab Nat. Phys. 6, 213 - 217 (2010). 
[15] Q. P. Unterreithmeier, E. M. Weig and J. P. Kotthaus, Nature 

458, 1001-1004 (2009). 
[16] A. Cleland, Nat. Phys. 5, 458-460 (2009). 
[17] A. Schliesser, O. Arcizet, R. Riviere, G. Anetsberger and T. J. 

Kippenberg, Nat. Phys. 5 509-514 (2009). 
[18] K. H. Lee, T. G. McRae, G. I. Harris, J. Knittel, and W. P 

Bowen Phys. Rev. Lett. 104, 123604 (2010) 
[19] P. F. Cohadon, A. Heidmann and M. Pinard, Phys. Rev. Lett. 83 

3174-3177 (1999); 
[20] M. Pinard, P. F. Cohadon, T. Briant and A. Heidmann, Phys. 

Rev. A 63 013808 (2000). 
[21] M. Hossein-Zadeh, and K. J. Vahala, Opt. Lett. 32 1611 (2007). 
[22] I. S. Grudinin, H. Lee, O. Painter, and K. J. Vahala, Phys. Rev. 

Lett. 104, 083901 (2010). 
[23] M. C. Cross, A. Zumdieck, R. Lifshitz and J. L. Rogers, Phys. 

Rev. Lett. 93 224101 (2004). 
[24] M. Spletzer, A. Raman, H. SumaU and J. P. Sullivan, Appl. 

Phys. Lett. 92 114102(2008); 
[25] M. Sato, B. E. Hubbard and A. J. Sievers, Rev. Mod. Phys. 78 

137-157 (2006). 
[26] R. J. Schoelkopf and S. M. Girvin, Nature 451 664-669 (2008). 
[27] D. Kleckner and D. Bouwmeester, Nature 444 75-78 (2006). 



[28] S. C. Masmanidis, R. B. Karabalin, L De Vlaminck, G. Borghs, 

M. R. Freeman and M. L. Roukes, Science 317, 780-783 

(2007). 
[29] K. L. Ekinci and M. L. Roukes, Rev. Sci. Instrum. 76, 061101 

(2005) 
[30] R. G. Knobel and A. N. Cleland, Nature 424 291-293 (2003). 
[31] N. E. Flowers-Jacobs, D. R. Schmidt and K. W. Lehnert, Phys. 

Rev. Lett. 98 096804 (2007). 
[32] M. D. LaHaye, O. Buu, B. Camarote, and K. C. Schwab, Sci- 
ence 304 74-77 (2004). 
[33] S. Etaki, M. Poot, I. Mahboob, K. Onomitsu, H. Yamaguchi and 

H. S. J. van der Zant, Nat. Phys. 4 785-788 (2008). 
[34] G. Binnig, C. R Quate and Ch. Gerber, Phys. Rev. Lett. 56, 

930-933 (1986). 
[35] U. Mohideen and A. Roy, Phys. Rev. Lett. 81 4549-4552 

(1998); 
[36] J. L. Hutter and J. Bechhoefer, J. Vac. Sci. Technol. B 81, 2251- 

2253 (1994); 
[37] D. Rugar, R. Budakian, H. J. Mamin & B. W. Chui, Nature 430 

329-332 (2004) 
[38] T. Briant, P F. Cohadon, A. Heidmann, and M. Pinard Phys. 

Rev. A 68, 033823 (2003). 
[39] O. Arcizet, C. MolineUi, T. Briant, P. F. Cohadon, A. Heid- 
mann, J. M. Mackowski, C. Michel, L. Pinard, O. Franais and 

L. Rousseau New J. Phys. 10, 125021 (2008). 
[40] D. K. Armani, T. J. Kippenberg, S. M. Spillane and K. J. Vahala, 

Nature 421 925-929 (2003). 
[41] C. W. Gardiner, and M. J. Collett, Phys. Rev. A 31, 3761-3774 

(1985) 
[42] Bowen, W. P. PhD Thesis, Australian National University, 

(2003). 
[43] T. Carmon, L. Yang and K. J. Vahala, Opt. Exp. 12 4742-4750 

(2004). 
[44] T. J. Kippenberg, H. Rokhsari, T. Carmon, A. Scherer, and K. 

J. Vahala, Phys. Rev. Lett. 95 033901 (2005). 
[45] A. Schliesser, G. Anetsberger, R. Riviere, O. Arcizet and T. J. 

Kippenberg, New J. Phys. 10 095015 (2008). 
[46] L. Landau and E. Lifshitz, Course of Theoretical Physics: Sta- 
tistical Physics Pergamon, New York, (1958), Chap. 12. 
[47] G. Anetsberger, R. Riviere, A. Schliesser, O. Arcizet and T. J. 

Kippenberg, Nat. Phot. 2 627 (2008). 
[48] O. Arcizet, R. Riviere, A. Schliesser, G. Anetsberger, and T. J. 

Kippenberg, Phys. Rev. A 80, 021803(R) (2009). 
[49] M. Li, W. H. P Pernice, C. Xiong, T. Baehr- Jones, M. Hochberg 

and H. X. Tang, Nature 456, 480-484 (2008). 
[50] T. Carmon, H. Rokhsari, L. Yang, T. J. Kippenberg, and K. J. 

Vahala, Phys. Rev. Lett. 94 223902 (2005). 



